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ABSTRACT 

The Naval Research Laboratory (NRL) Layered Ocean Model (NLOM) with an embedded bulk-type mixed 
layer model is used to examine the effects of ocean turbidity on sea surface temperature (SST) and ocean mixed 
layer depth (MLD) simulations over the global ocean. The model accounts for ocean turbidity through depth- 
dependent attenuation of solar radiation in the mixed layer formulation as determined from the diffusive atten¬ 
uation coefficient at 490 nm (k^^) obtained by the Sea-Viewing Wide Field-of-View Sensor (SeaMFS). Inter- 
annual model simulations are used to assess the first-order effects of ocean turbidity on SST and MLD simulation. 
Results are reported from three model experiments performed using different values for the attenuation of 
photosynthetically available radiation (^par)- It is shown that, although allowing incoming solar radiation to vary 
in time and space is desirable for predicting SST, in an OGCM use of a constant ApAp with a value of 0.06 m"‘ 
is generally sufficient in the deep ocean. The daily averaged SST time series from the three NLOM simulations 
are verified against daily in situ SSTs reported from 12 moored buoys in 1996 and 1997. Model results show 
that allowing the possibility of solar heating below the mixed layer reduces the root-mean-square error (rmse) 
difference between the daily yearlong model and buoy SST time series by up to 0.4°C and r^uces the rmse at 
11 of the 12 buoy locations. Although using spatially and temporally varying fcpAp versus a constant Atpar - 0.06 
m"* (which is representative over most of the globid ocean) had low impact overall, using it generally reduced 
the rmse at low latitudes, and using it can have a substantial impact locdly in space and time. The model MLD 
results show low sensitivity to the ^par value used. 


1. Introduction 

Ocean mixed layer models play an important role in 
simulated air-sea interactions and are required to ex¬ 
plain upper-ocean characteristics (e.g., Martin 1985; 
Sterl and Kattenberg 1994; Schopf and Loughe 1995). 
The major difficulty in parameterizing the ocean mixed 
layer is that the ocean boundary layer is not as fully 
observed as the atmospheric boundary layer (e.g.. Smith 
et al. 1996). Thus, specific attention must be given to 
the atmospheric forcing that drives the upper ocean. 

Sensible and latent heat fluxes, longwave and short¬ 
wave radiation, wind stress, and ocean surface temper- 
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ature constitute the principal means for coupling be¬ 
tween the atmosphere and ocean. Among these vari¬ 
ables, the vertical distribution of shortwave radiation 
within the ocean is significant for upper-ocean studies 
because radiation can penetrate below the mixed layer 
and destabilize the stratification (e.g., Austin and Pet- 
zold 1986; Morel and Antonie 1994; Morel and Mari- 
tonera 2001), thereby affecting SST during upwelling/ 
downwelling events (Nakamoto et al. 2000). Knowledge 
of the shortwave radiation within the ocean is particu¬ 
larly critical for understanding dynamical and thermo¬ 
dynamical processes occurring in the ocean surface 
mixed layer. For example, the amount of solar irradiance 
affects sea surface temperature (SST) and ocean surface 
mixed layer depth (MLD) by changing the surface en¬ 
ergy budget of the water (e.g., Gallimore and Houghton 
1987; Lewis et al. 1990), and it controls the rates of 
photosynthetic processes (e.g., Frouin et al. 1989; Brock 
and McClain 1992; Brock et al. 1993). We note here 
that unlike the atmospheric boundary layer, which is 
nearly transparent, the ocean boundary layer strongly 
absorbs solar radiation. Thus, changes in upper-ocean 
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characteristics are closely tied to variations in solar ra¬ 
diation within the ocean and its attenuation within the 
water column, which need to be taken into account in 
numerical ocean modeling studies. Radiation attenua¬ 
tion might also be significant for long-term climate var¬ 
iability if the attenuation depth scale is greater than the 
MLD. For example, a fraction of the solar radiation 
heating might be absorbed below the thermocline during 
the onset of an El Nino event (e.g., Woods 1994). This 
heating may not be redistributed easily by the mixed 
layer and could be significant on interannual time scales. 

Bulk mixed layer models (e.g., Garwood 1977; Niller 
and Kraus 1977; Price et al. 1986) and diffusion-based 
mixed layer models (e.g., Mellor and Yamada 1982) 
have different responses to solar heating as reviewed by 
Kantha and Clay son (1994). The sensitivity to the dis¬ 
tribution of the solar radiation (Simpson and Dickey 
1981; Dickey 1983; Le Treut et al. 1985) is primarily 
because vertical stratification plays a controlling role in 
the vertical turbulent mixing. As was shown in Simonot 
and Le Treut (1986), part of the solar radiation can be 
absorbed below the MLD at low latitudes. Therefore, 
the modification of the water optical properties can 
change the density field of the ocean and the dynamical 
response of the oceanic layers to surface wind stress 
forcing (e.g., Kara et al, 2000a). High turbidity also 
raises the temperature in a shallow mixed layer, thus 
changing the SST and stratification. This implies that a 
good representation of optical properties of seawater is 
necessary for accurate modeling of SST and MLD. The 
importance of solar extinction in the upper ocean in 
relation to the diurnal cycle was discussed using a mixed 
layer model based on second-moment closure of tur¬ 
bulence (Kantha and Clayson 1994). They explained 
that the rise in SST for a very shallow MLD is highly 
sensitive to the parameterization of the solar extinction 
that includes different Jerlov (1976, 1977) water types. 
When using a clear water type instead of a variable water 
type in their simulations, the diurnal cycle of SST es¬ 
sentially disappeared. Ofiier studies (e.g.. Woods et al. 
1984; Lewis et al. 1990) have also discussed the effects 
of optical water types on the ocean surface mixed layer. 
They indicated the need to consider the penetration of 
solar irradiance in developing simulations of tropical 
Pacific SST. 

The sensitivity of SST and MLD evolution to solar 
transmission in large-scale ocean general circulation 
model (OGCM) studies has also been investigated. For 
example, Schneider and Zhu (1998) showed with a cou¬ 
pled atmosphere-ocean GCM that inclusion of sunlight 
penetration was essential to correctly simulate the strong 
semiannual SST cycle in the eastern equatorial Pacific. 
This was achieved using a constant penetration depth 
of 15 m and did not use any satellite-based ocean color 
data to account for water turbidity. Murtugudde et al. 
( 2002 ) illustrated the importance of a realistic treatment 
of penetrating radiation in climate models by using spa¬ 
tially varying annual-mean attenuation depths that were 


estimated from remotely sensed chlorophyll concentra¬ 
tions provided by the Coastal Zone Color Scanner 
(CZCS). While their OGCM study was limited to the 
Tropics because of the available coverage of CZCS, the 
study did show significant improvement in model SST 
simulations when using annual-mean attenuation depths 
versus a constant value of 17 m. The study by Nakamoto 
et al. ( 2001 ) examined the effects of including space- 
and time-varying solar attenuation with depth by using 
monthly CZCS chlorophyll concentrations where avail¬ 
able over the World Ocean. Their analysis of the model 
simulations focused on the equatorial Pacific. Ohlmann 
et al. (1996) used monthly CZCS data with a fully spec¬ 
tral radiation model to show that solar radiation pene¬ 
tration can be a significant term in the mixed layer heat 
budget for the tropical regions. Rochford et al. (2001) 
constructed a global monthly mean climatology for the 
attenuation of photosynthetically available radiation 
(PAR), denoted as ^p^R, using a diffusive attenuation 
coefficient at 490 nm (^ 490 ). They showed that the in¬ 
clusion of solar subsurface heating is important for im¬ 
proving the global OGCM simulation of SST in the 
equatorial regions, while subsurface heating is not im¬ 
portant for MLD simulations. Although these studies 
indicate the importance of including a penetrating ra¬ 
diation scheme to improve model SST and MLD sim¬ 
ulations, none examined the effects of water turbidity 
on SST and MLD simulations on short time scales (e.g., 
daily). 

The main goal of this paper is to indicate the regions 
where water turbidity variability can affect SST and 
MLD simulations over the global ocean on daily time 
scales. For this purpose we examine daily SST and MLD 
sensitivity to subsurface heating using interannual sim¬ 
ulations from the Naval Research Laboratory (NRL) 
Layered Ocean Model (NLOM) with an embedded 
mixed layer. All model simulations presented in this 
paper use atmospheric forcing with high temporal (6 
hourly) resolution. Model-data comparisons are per¬ 
formed using daily moored buoy observations. 

This paper is organized as follows. In section 2, a 
brief description is given of the NLOM mixed layer 
implementation, surface forcing, turbulence model, and 
model parameterizations. In section 3, the effects of 
ocean turbidity on the model simulations are explained 
by comparing the daily model SST time series with daily 
buoy SST time series over the different regions of the 
global ocean. Some comparisons are also made between 
the findings obtained from NLOM and the findings of 
similar studies that used different OGCMs. A summary 
and conclusions of this paper are presented in section 4. 

2. Mixed layer model description 

The numerical ocean model used is a primitive equa¬ 
tion layered formulation in which the equations have 
been vertically integrated through each layer. The 
NLOM is a descendent of the model by Hurlburt and 
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Fig, 1. Schematic illustration of mixed layer structure in the NLOM 
adapted from Wallcraft et al. (2003). For simplicity, only three ther¬ 
modynamic layers and the mixed layer are shown, although six layers 
plus the mixed layer are used in this study. Prognostic variables are 
layer density (p^X layer thickness (h^X layer transport per unit width 
(h^y^X mixed layer temperature (r„), and MLD {h^X Atmospheric 
forcing is wind stress (t), air temperature (ro), mixing ratio at 10 m 
above the sea surface {qj, longwave radiation (Glw)» and downward 
solar irradiance at the sea surface (Gsol)- 


Thompson (1980) with enhancements by Wallcraft 
(1991), Wallcraft and Moore (1997), and Moore and 
Wallcraft (1998). It can be run in reduced-gravity mode, 
in which the lowest layer is infinitely deep and at rest, 
or in finite-depth mode, which allows realistic bottom 
topography. In addition, NLOM can be run in hydro- 
dynamic mode (spatially and temporally constant den¬ 
sity within each layer; e.g., Hurlburt et al. 1996) or in 
thermodynamic mode (spatially and temporally varying 
density within each layer—i.e., density is a prognostic 
variable; e.g,, Metzger and Hurlburt 2001). The model 
boundary conditions are kinematic and no slip. Prog¬ 
nostic variables are layer density, layer thickness, and 
layer volume transport. The model includes entrainment 
(detrainment) processes by allowing the exchange of 
mass, momentum, and heat between the layers. 

The model domain used for this study covers the 
global ocean from 72®S to 65°N, gridded to a resolution 
of 0,5° in latitude and 0.703 125° in longitude. The 
lateral boundaries follow the 200-m isobath. One of the 
major advantages of NLOM over other types of OGCMs 
such as z-level and sigma-coordinate models is its lower 
computational cost for the same model domain and hor¬ 
izontal resolution. One reason is that we can use lower 
vertical resolution to realistically represent the ocean 
circulation. For example, in the (W) model there are 


only six dynamical layers plus a semipassive mixed lay¬ 
er in the vertical with layer thickness that varies in both 
time and space. NLOM is also a single efficient portable 
and scalable computer code that can run any of the 
model configurations on a variety of computing plat¬ 
forms (Wallcraft and Moore 1997). 

In this paper, we give only a brief description of the 
embedded mixed layer parameterizations. See Wallcraft 
et al. (2003) for a full description of the model for¬ 
mulations and the W global model used here. See Kara 
et al. (2000b, 2002a) for the parameterizations used in 
calculating atmospheric forcing fields. 


a. Mixed layer equations 

The NLOM has typically been run without an explicit 
mixed layer, which is equivalent to assuming the mixed 
layer is always inside the upper layer, but now the model 
has been extended with an “almost passive” embedded 
well-mixed surface turbulent boundary layer. A sche¬ 
matic illustration of the global NLOM with an embed¬ 
ded mixed layer is shown in Fig. 1. 

The embedded mixed layer model employed here car¬ 
ries prognostic equations for the SST (T^) and MLD 
(/i„,) as follows; 

SST: 


^ + V. • vr„ = - t ,) 


MLD: 


^ Q p-hjhp X 

PoLpartm "m 

( 1 ) 


d(hj 

dt 




( 2 ) 


The symbols appearing in the equations are defined in 
the appendix along with their units. Major free param¬ 
eters in Eqs. (1) and (2) are net surface heat flux (QJ, 
temperature shear at the base of the mixed layer [AJ^ 
= (T„ - ATJ- rj, and vertical mixing velocity ((oj. 
These free parameters are obtained from a surface en¬ 
ergy budget (section 2b) and from a continuous model 
temperature profile, respectively. The temperature 
change across the mixed layer is Ar,„, specified as a 
function of latitude based on the NMLD monthly cli¬ 
matology (Kara et al. 2002b). The values of A7„ range 
from 0.1°C! at high latitudes to 1.5°C at low latitudes. 

In the model simulations, the Laplacian temperature 
diffusion in Eq. (1) is typically turned off (i.e., = 

0) because the van Leer monotonic scheme (Lin et al. 
1994) used for advection contributes sufficient nonlinear 
diffusion for stability. This scheme allows NLOM to be 
run using a very small minimum imposed MLD (e.g., 
10 m). 
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b. Surface energy balance 

The net surface heat flux that has been absorbed (or 
lost) by the upper ocean to depth z, Q(zX is parameter¬ 
ized as the sum of the downward surface solar irradiance 
{2sol)» upward longwave radiation (Qlw)* and the 
downward latent and sensible heat fluxes (2^ and Qs, 
respectively). The solar irradiance at the ocean surface 
ranges in wavelength from about 300 to 2800 nm and 
is composed of three general regions: the ultraviolet 
(UV) below 400 nm, the visible at 400-700 nm, and 
the infrared (JR) above 700 nm. The light available for 
photosynthesis by photoplankton is the PAR and is de¬ 
fined as the 350-700-nm range of the spectrum (e.g., 
Liu et al. 1994). It accounts for 43%-50% of the solar 
irradiance at the sea surface (Rochford et al. 2001). 

For the NLOM surface forcing, the surface solar ir¬ 
radiance is decomposed into its infrared ((2 ir) arid Pen¬ 
etrating radiation (Qp) as Qsol = Qm + Qpy thereby 
enabling us to write 

Q(z) = jSiR + Qpiz) ■“ 0LW Ql Qs* (^) 

Since the component is absorbed within the first few 
centimeters of the sea, a depth which is much less than 
the minimum MLD imposed in most OGCMs (10 m for 
the NLOM), all of the IR radiation arriving at the air- 
sea interface is used to heat the surface layer of the 
ocean model. The Qp component for the PAR and UV 
portion of the spectrum that penetrates to greater depths 
can be well represented by a single exponential. Based 
on global monthly climatologies of Qp and Qsol^ Qir “ 
0.49gsoL is assumed everywhere over the global ocean. 
Thus, the expression for subsurface heating is obtained 
as follows: 

Q(z) = (4) 

hp 1/^PAR» 3nd (5) 

Qa “ GsOL "" QlW Qs^ (6) 

Here Q^ is the net heat flux at the ocean surface. The 
rate of surface heating/cooling of the mixed layer is 
simply obtained by evaluating this expression at the 
MLD (i.e., z = hj), and the remaining solar radiation, 
Qih^X is applied below the mixed layer. The radiation 
absorption length scale in Eq. (5) is expressed as hp = 

It is noted that we use Eq. (4) although fully spectral 
representations are available (e.g.. Morel and Marito- 
rena 2001). In general, OGCMs that have few layers 
near the surface (e.g., Schneider and Zhu 1998), with 
a minimum MLD, for example, 10 m, use simple solar 
irradiance approximations, such as single or bimodal 
exponential parameterizations (Paulson and Simpson 
1977; Zaneveld and Spinrad 1980). This is because, for 
depths greater than 10 m, the penetrative solar flux can 
be accurately determined by resolving just the 300-745- 
nm spectral region, which is well represented by a single 
exponential below 10 m. 


As explained in detail in Kara et al. (2003, unpub¬ 
lished manuscript, hereinafter KHR), the attenuation 
depth for PAR differs depending on water type (Simonot 
and Le Trent 1986) and is quite variable over the global 
ocean. The extent to which PAR penetrates into and 
below the mixed layer has been quantified using re¬ 
motely sensed satellite color data and used in other 
OGCM studies (e.g., Schneider and Zhu 1998; Murtu- 
gudde et al. 2002; Kara et al. 2003b). To properly in¬ 
clude effects of turbidity, the model therefore reads in 
monthly fields (KHR) based on Sea-Viewing Wide 
Field-of-View Sensor (SeaWiFS) data (McClain et al. 
1998). Thus, our approach amounts to a one-band versus 
two-band approach, where the red and near-infrared ra¬ 
diation is completely absorbed within the minimum 
MLD of 10 m imposed by the NLOM, and the pene¬ 
trating radiation portion is allowed to have space and 
time variations. 

Latent and sensible heat fluxes at the air-sea interface 
are calculated using efficient and computationally in¬ 
expensive simple bulk formulas that include the effects 
of dynamic stability (Kara et al. 2002a). The combi¬ 
nation of accuracy and ease of computation of this meth¬ 
od makes it the one preferred for computing air-sea 
fluxes in the NLOM. Note that both sensible and latent 
heat fluxes are calculated using mixed layer temperature 
(T„) at the model time step. Radiation flux (shortwave 
and longwave fluxes) is so dependent on cloudiness that 
this is taken directly from European Centre for Medium- 
Range Weather Forecasts (ECMWF) data (ECMWF 
1995; Gibson et al. 1997) for use in the model. Basing 
fluxes on the model SST automatically provides a phys¬ 
ically realistic tendency toward the “correct” SST. If 
the model SST is too high or low, the flux is reduced 
or increased relative to that from the correct SST. In 
NLOM, the accuracy of the SST is also enhanced by 
relaxing the dynamic layer densities below the mixed 
layer back toward climatology (monthly in layer 1, an¬ 
nually otherwise), in addition to applying the heat flux 
(Kara et al. 2003b). There is no direct SST relaxation 
term in the equation, but entrainment at the base of the 
mixed layer allows the dynamical layer density relax¬ 
ation to influence SST. 

c. Turbulence model 

The rate of mixed layer deepening or retreat, (o,„, is 
determined using a modified version of the KT model 
(Kraus and Tbmer 1967; Niiler and Kraus 1977). This 
involves solving the vertically integrated turbulent ki¬ 
netic energy (TKE) equation for a stationary budget. 
The TKE budget (P) implemented in the NLOM for the 
present study is already derived in KHR with particular 
attention to the attenuation of solar radiation with depth: 
P = (W3 - mi)ul 

- - m^Sb^w^ + mju^ + (7) 

where 
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ul 




8 

f 


MA 

Po ’ 

(8) 

- Qp), 

Po^pw 

(9) 

1, for b^w^ < 0 

(10) 

0, for S: 0, 

inax(|2fl sm(^)|, /+), 

(11) 


and contains the contribution due to attenuation of 
solar radiation that will be specified further below. 

When shallowing (i.e., P < 0), mixing occurs toward 
the equilibrium depth as follows: 


h* = 


(OTj - mi)u% 


B ^ 

- rnf,8b^w^ + mju% 

’ «« 


and (12) 


(o,^ = K) “ hj, (13) 

A relaxation time scale (cr^) is introduced to account 
for a delayed retreat of the mixed layer to the equilib¬ 
rium depth. When deepening (i.e., P > 0), the available 
TKE is converted to potential energy as follows: 

P 

” ga{TJ max(Ar„ 

AT, = (r„ - Arj - T,. (15) 

A minimum value is imposed on the MLD {h^ = 10 
m) and AT* represents the temperature difference at the 
base of the mixed layer. Since A 7^ is not guaranteed to 
be positive, an alternative minimum value (AT^ = 
0.2°C) is specified. For further information on the details 
and rationale for the above the reader is referred to 
Wallcraft et al. (2003). 

Various possible forms can be chosen for the contri¬ 
bution due to radiation penetration (B„d) and these can 
have functional dependencies upon that are highly 
nontrivial. To retain a reasonable level of simplicity so 
that insight can be gained when including and excluding 
the effects of subsurface heating, we employ the fol¬ 
lowing for Brad ill this study: 


(16) 

2 PoSw 

This choice reproduces the essential result for when 
hp becomes extremely small, 

Hm (17) 

while capturing the maximum possible effect of the sub¬ 
surface heating for finite values of hp. As noted in KHR, 
the form derived for B^ad is strictly valid only within the 
surface layer, which does not necessarily extend over 
the entire depth of the mixed layer, and serves only to 
guide the choice of parameterization. With this chosen 


form for B^^d one has the property that B^^d ^ 0 {Qp ^ 
0), thereby facilitating interpretation of the effects of 
including subsurface heating. For example, when shal¬ 
lowing (i.e., P < 0), the equilibrium depth now becomes 


hi = 


(m3 - m^)ul 


(18) 


“ m^db^w^ + mjul 


where 


= ^^(2. - Gpc-*"'*-). (19) 

Po^pw 


For a finite value of hp, the buoyancy flux b^u^ will be 
decreased relative to that when hp is infinitesimally 
small (equivalent to complete absorption of the heat flux 
at the air-sea interface), thereby implying a deeper 
mixed layer when subsurface heating is included in an 
OGCM with an embedded mixed layer. Similarly, when 
deepening (i.e., P > 0), the mixed layer will be deeper 
with the inclusion of subsurface heating (finite vs in¬ 
finitesimal hp) because b^u^ will become smaller and 
increase the amount of TKE (B) available for entraining 
fluid into the mixed layer. 


d. Model spinup 

The six-layered thermodynamic model without the 
mixed layer is spun up to statistical energy equilibrium 
and then continued for 5 years with a mixed layer and 
climatological forcing. Climatological monthly means 
of the thermal forcing are obtained from the Compre¬ 
hensive Ocean-Atmosphere Data Set (COADS) (da Sil¬ 
va et al. 1994). Thermal forcing includes shortwave (in¬ 
coming solar) plus longwave radiation {Qp\ air tem¬ 
perature (7^) at 10 m, and the air mixing ratio {q^) at 
10 m. Scalar wind speed (v„) is obtained from the input 
wind stress, which has 6-hourly variability. All clima¬ 
tological model simulations with a mixed layer are per¬ 
formed using climatological 6-hourly hybrid winds. 
These consist of monthly Hellerman and Rosenstein 
(1983) wind stresses (HR) plus ECMWF wind anom¬ 
alies (Wallcraft et al. 2003). We add a high-frequency 
component to the climatological forcing because of its 
impact on the mixed layer and because the model is 
targeted at simulations forced by high-frequency inter¬ 
annual atmospheric fields from operational weather cen¬ 
ters. After reaching statistical energy equilibrium, the 
climatological model simulation with a mixed layer is 
extended using interannual surface forcing from 1996 
and 1997. In this case, the 6-hourly thermal forcing from 
ECMWF includes the shortwave plus longwave radia¬ 
tion {Qp), air temperature (7^) at 10 m, and air mixing 
ratio {q^) at 10 m. None of the model simulations pre¬ 
sented in this paper include assimilation of any SST 
data. 
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Fig. 2. Percentage absorption of solar radiation with depth (h) 
calculated using a single exponential decay formula. These profiles 
are good representatives of solar transmission in the visible band 
within the mixed layer below 10 m. The same approach is also used 
in NLOM. 


3. Turbidity effects on NLOM simulations 

In this section, sensitivity of model results to water 
turbidity is examined with a particular focus on SST 
and MLD simulations. The model results are then com¬ 
pared with daily SST and MLD time series obtained 
from buoys located in different regions of the global 
ocean. 

a. Global turbidity ejfects 

Before analyzing the model results, we first examine 
where turbidity might be important to the annual mean 
and seasonal cycle of the mixed layer. Our purpose is 
to indicate the regions where variability can affect 
the SST simulation over the global ocean. The greatest 
interest is in those situations where subsurface heating 
can occur below the mixed layer, thereby resulting in a 
cooler SST than is obtained by assuming complete ab¬ 
sorption of solar irradiance within the bulk mixed layer. 
This is expected to occur in regions of small kp^R values. 
The reverse situation of large values is expected to 
produce complete (or almost complete) absorption of 
the solar radiation within the mixed layer because a 
minimum MLD of 10 m is imposed by the NLOM. It 
is noted that, based on a single exponential approach 
[100-49.0 exp(-/ikpAR)] as shown in Fig. 2, all large 
itpAR values (i.e., kp^R > 0.20 m"0 are essentially uni¬ 
form. 

Three different kp^p values are considered: 1) a 
monthly itp^R dataset to account for the seasonal and 


spatial variation of subsurface heating, 2) a global con¬ 
stant value of itpAR = 0.06 m"^ that is representative of 
most open-ocean conditions, and 3) an unrealistically 
large value of ^p^R = 99 m‘^ to ensure the complete 
absorption of solar radiation within the mixed layer. 
Case 2 corresponds to an ^-folding penetration depth of 
hp = itpAR ” 16.7 m with 95% absorption by 50-m depth. 
This is close to hp values used in other studies: 15 
(Schneider and Zhu 1998), 17 (Murtugudde et al. 2002), 
and 23 m (Nakamoto et al. 2001). Case 3 represents 
traditional OGCM approaches (e.g., Yuen et al. 1992; 
Murtugudde et al. 1995). 

For all three cases the monthly mean radiation ab¬ 
sorbed below the mixed layer is calculated using Qp 
exp(“z//ip), with hp [cf. Eq. (5)] defined as a temporally 
interpolated monthly global field. These calculations use 
global climatologies of monthly ^p^R and monthly op¬ 
timal MLD based on observations (Kara et al. 2003a). 
Then, the differences of cases 1-2 and cases 1-3 
are calculated month by month and presented as a per¬ 
centage relative to the net surface heat flux at the ocean 
surface Note here that case 1 - case 2 represents 
the percentage difference absorbed from using monthly 
versus kp^ = 0.06 m“*, and case 1 - case 3 is just case 
1 because Qp is assumed to be zero for case 3. Last, 
global maps of the annual mean and the maximum value 
over the 12 months are formed. The reason for using a 
worst-month approach is the root-mean-square error 
(rmse) would not be a useful measure of turbidity impact 
because the error is likely to be near zero whenever the 
MLD is deep. The annual-mean climatological MLD is 
as shallow as 10 m in some regions where absorption 
is greater, such as parts of the equatorial region and the 
North Pacific (see KHR). This is especially evident in 
summer months. For example, strong solar heating oc¬ 
curs in conjunction with optical depths on the order of 
10-15 m throughout the year in the equatorial Pacific. 

The impact of a space- and time-varying 4 ar can be 
most easily determined from the ratio of PAR below the 
MLD relative to the PAR at the air-sea interface. An 
examination of the annual mean difference between ex¬ 
periments for this PAR fraction expressed as a per¬ 
centage (Figs. 3a and 3b) reveals that ^p^R has essentially 
no effect over most of the global ocean (regions dis¬ 
played as white). This occurs where PAR penetration 
below the MLD is low because of a relatively deep 
mixed layer most of the time. The relatively deep annual 
mean MLD, as shown in KHR, tends to correlate well 
with Fig. 3a and/or high turbidity, especially in relation 
to Fig. 3c. Using a monthly kp^R makes a large difference 
(in at least one month) at high latitudes and in the equa¬ 
torial ocean (Fig. 3c). The sensitivity to solar attenuation 
in the Tropics agrees with earlier studies (Murtugudde 
et al. 2002; Nakamoto et al. 2001; Schneider and Zhu 
1998). Note that the absorption anomaly only indicates 
when solar heating within the mixed layer is signifi¬ 
cantly affected by turbidity and does not indicate its 
influence on SST. 
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Fig. 3. Changes in PAR below the MLD as percentage of PAR at the surface: (a) annual-mean difference when using a monthly varying 
JtpAR relative to complete PAR absorption within the mixed layer (i.e., the effect of ignoring ^par), (b) annual-mean difference when using 
a monthly varying JfepAR relative to itpAR = 0,06 m-^ and (c) largest difference of all 12 months for case (b), indicating where using a monthly 
itpAR can make a large difference in model simulations. 


Based on the observational results presented above, 
interannual model simulations that use three different 
jfcpAR values are performed to investigate the effects of 
ocean turbidity for the model SST and MLD simulations 
in 1996 and 1997 (see Table 1 for a brief explanation 
of each experiment performed). The model spinup and 
forcing were already described in section 3. For exper¬ 


iment 1, spatially and monthly varying values in¬ 
terpolated to the global NLOM grid are used. For ex¬ 
periment 2, the global ocean turbidity is set to a constant, 
^AR — 0-06 m“L In experiment 3, all of the solar ra- 
iation is absorbed in the mixed layer by using a very 
large ^par value. Differences between experiments 1 and 
2 are only expected near the equator and in coastal 
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Table 1. The three interannual NLOM simulations for 1996 and 
1997. The model has an embedded bulk-type mixed layer. It used 
wind and thermal forcings (i.e., air temperature at 10 m, air mixing 
ratio at 10 m, shortwave and longwave radiation) from the 6-hourIy 
ECMWF data. 


Expt 

^PAR 

Description of the expt 

i 

Variable 

Spatial and temporal turbidity 

2 

0.06 m-‘ 

Constant turbidity over the global ocean 

3 

99 m-‘ 

All solar radiation absorbed in mixed layer 


regions where larger space- and time-varying val¬ 
ues occur (Fig. 3). Differences in SST between these 
two simulations only occur in regions where a shallow 
MLD coincides with a monthly varying 4 ar that differs 
substantially from 0.06 m"^ 


b. NLOM model-data comparisons 

For model-data comparisons, the daily averaged SSTs 
obtained from the three interannual NLOM simulations 
are compared with daily averaged SST time series from 
moored buoys located in different regions of the global 
ocean. These are obtained from the Tropical Atmosphere- 
Ocean (TAO) array (McPhaden 1995) and the National 
Oceanic Data Center (NODC). In both datasets, SST is 
measured at a depth of 1 m below the sea surface. The 
model and buoy SSTs are compared using several sta¬ 
tistical metrics (e.g., Stewart 1990; Murphy 1988) as de¬ 
scribed next. Let v, (f = 1,2,..., N) be the set of N 
observed (i.e., buoy) SSTs, and let e, (f = 1,2,..., N) 
be the set of NLOM SSTs. Also let v(J) and a’„(a;) be 
the mean and standard deviations of the buoy (NLOM) 
SSTs. The statistical measures can then be expressed as 
follows: 


ME — e — V, 

^ , _-^(^/ “ ^) 

R = (Vi- v)- - 

n Si (o-^CTe) 

1 ” 

(rmse)2 = - 2 (e, v^y, and 
n 


SS = 1 - 


(vmsef 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 


where ME is the mean error, rmse is the root-mean- 
square difference, R is the correlation coefficient, and 
SS is the skill score. Because the standard deviation of 
SST is different for each buoy, rmse between the model 
and the buoy will not entirely reflect model perfor¬ 
mance. Thus, we also include a dimensionless measure 
(SS) that takes account of the bias between the model 
and buoy SST; SS = 1 for a perfect estimate (Murphy 
1988). Negative SS values indicate the model SSTs do 
not compare well with the buoy SSTs. 

Table 2 statistically compares the daily averaged 
NLOM SSTs from the three experiments with TAO buoy 
SSTs in 1997, and Table 3 shows the same for NODC 
buoy SSTs outside the equatorial region. Overall, the 
ME values are small between the NLOM and buoy SSTs 
for all locations. The NLOM is able to capture SST 
variability well because correlation coefficients are close 
to /? = 1 in most of the cases. Large and positive SS 
values show that the NLOM is able to simulate SST 
with skill. This is especially true when turbidity effects 
(i.e., variable ^p^r) are taken into account (expt 1). Using 
a constant fcp^R (expt 2) rather than variable did not 
significantly affect the results for the locations used in 
this study. However, we clearly see that absorption of 
all solar radiation into the mixed layer (expt 3) simulates 
the SST less accurately. This is not surprising because 


Table 2. Statistical verification of daily SST time series between NLOM and TAO buoys at different locations in the equatorial ocean. 
Comparisons are made during a 1-yr period (1997). All statistical calculations are based on 365 daily SST values. 


Location 

Expt 

Rmse CC) 

ME CQ 

^BUOY ( O 

^NLOM ( 

R 

SS 

(0°, 110°W) 

1 

1.31 

-0.26 

1.45 

0.99 

0.50 

0.18 

2 

1.35 

-0.38 

1.45 

1.02 

0.50 

0.14 


3 

1.61 

0.07 

1.45 

1.14 

0.34 

-0.23 

(0^ 125°W) 

1 

1.21 

0.05 

1.64 

1.06 

0.68 

0.46 

2 

1.23 

0.12 

1.64 

1.09 

0.67 

0.44 


3 

1.35 

0.31 

1.64 

1.20 

0.61 

0.33 

(0°, 155°W) 

1 

0.77 

-0.13 

1.44 

1.03 

0.86 

0.72 

2 

0.78 

-0.13 

1.44 

0.97 

0.87 

0.71 


3 

0.87 

0.09 

1.44 

1.09 

0.80 

0.64 

(0°, 170°W) 

1 

0.62 

-0.32 

0.89 

0.72 

0.80 

0.51 

2 

0.66 

-0.40 

0.89 

0.70 

0.80 

0.44 


3 

0.64 

-0.06 

0.89 

0.89 

0.75 

0.49 

(5°N, 155°W) 

1 

0.55 

-0.14 

0.98 

0.97 

0.85 

0.69 

2 

0.53 

-0.15 

0.98 

1.00 

0.87 

0.70 


3 

0.61 

-0.06 

0.98 

1.05 

0.82 

0.62 

(8'S, 110°W) 

1 

0.70 

-0.43 

0.71 

0.71 

0.68 

0.08 

2 

0.70 

-0.44 

0.71 

0.70 

0.68 

0.08 


3 

0.77 

-0.31 

0.71 

0.78 

0.62 

-0.17 
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Table 3. As in Table 2 but for 1996 and with statistical verification performed using NODC buoys outside the equatorial ocean. 


Location 

Expt 

Rmse (°C) 

ME CC) 

^^BUOY ( f') 

^^NLOM ( 

R 

SS 

(26X 94°W), Gulf of 

1 

0.58 

-0.36 

3.19 

3.21 

0.99 

0.97 

Mexico 

2 

0.61 

-0.40 

3.19 

3.19 

0.99 

0.96 

3 

0.63 

-0.01 

3.19 

3.61 

0.96 

0.96 

(29®N, 79°W), Florida coast 

1 

0.83 

0.41 

2.15 

2.49 

0.96 

0.85 

2 

0.80 

0.38 

2.15 

2.47 

0.96 

0.86 


3 

1.19 

0.69 

2.15 

2.79 

0.95 

0.70 

(38®N, 130°W), off southwest 

1 

0.97 

-0.42 

1.63 

2.07 

0.91 

0.64 

United States 

2 

0.97 

-0.49 

1.63 

2.05 

0.92 

0.64 

3 

1.03 

-0.28 

1.63 

2.18 

0.91 

0.60 

(41°N, 137°W), off northwest 

1 

0.67 

0.08 

2.94 

2.99 

0.97 

0.95 

United States 

2 

0.63 

.0.03 

2.94 

2.97 

0.98 

0.95 


3 

0.99 

0.32 

2.94 

3.29 

0.96 

0.89 

(46°N, 131°W), off northwest 

1 

0.83 

-0.01 

2.85 

3.09 

0.96 

0.91 

United States 

2 

0.77 

-0.08 

2.85 

3.04 

0.97 

0.93 

3 

1.03 

0.15 

2.85 

3.29 

0.96 

0.87 

(52°N, 156°W), Alaska coast 

1 

0.78 

0.41 

2.66 

2.87 

0.97 

0.91 

2 

0.67 

0.33 

2.66 

2.77 

0.98 

0.94 


3 

0.82 

0.43 

2.66 

2.89 

0.97 

0.91 


this choice can affect SST in two ways: 1) by allowing 
too much solar radiation to be absorbed in the mixed 
layer and 2) by influencing the stability of the upper 
ocean and thus SST. 


To further examine the effects of turbidity the ob¬ 
served SSTs from two moorings are compared with 
those obtained from the three NLOM experiments (Fig. 
4): one located in the central equatorial Pacific (0°, 



Fig. 4. Daily averaged SST time series comparisons between two buoys and three NLOM 
experiments (see Table 1). Two different buoy locations are used in two different years: (a) 
equatorial ocean (0°, i55°W) in 1997 and (b) off the northwest coast of the United States (46°N, 
131*^) in 1996. See Fig. 5 for monthly values at these locations. 
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Fig. 5. Climatological monthly mean of kp^R at the locations given 
in Tables 2 and 3. Note that jc-axis labels are written for every other 
month. 


155°W), and the other off the U.S. northwest coast 
(46®N, 13rW). Recall that experiment 1 is performed 
using turbidity at each buoy location as represented by 
the monthly SeaMFS 4 ar values. For the buoy at (0°, 
155°W), it is obvious that complete absorption of solar 
radiation within the mixed layer (expt 3) produces the 
highest SST from February through April and from Oc¬ 
tober through November. Note that 4 ar values at this 
location are greater than 0.06 m"* for February-April 
but not October-November (Fig. 5). Clearly, absorption 
of all solar radiation (expt 3) in the mixed layer causes 
a relatively large warm bias between the NLOM SST 
and buoy SST at this location. A similar warm bias was 
found in the annual mean SST by Schneider and Zhu 
(1998) between simulations assuming no solar penetra¬ 
tion and a constant attenuation depth for solar penetra¬ 
tion. As for (46°N, ISTW), the SST differences between 
experiment 1 (the standard simulation) and experiment 
2 are again small. These results agree with other studies 
(e.g., Murtugudde et al. 2002) that find larger SST 
changes occurring at low latitudes. Experiment 3 pro¬ 
duces the highest SST, and the difference with respect 
to experiments 1 and 2 is largest for June-September. 

To determine how the differences between the SST 
simulations may be due to differences in MLD, we used 
daily averaged subsurface temperature measurements 
from (0°, 155°W) to determine the daily MLD. Sub¬ 
surface parameters from the buoys include daily water 


temperatures at 10 discrete depths between 25 and 500 
m. For this buoy they are 25, 50, 75, 100, 125, 150, 
200, 250, 300, and 500 m. Subsurface temperature ac¬ 
curacy is about 0.01®-0.05°C for the different sensors 
used (Freitag et al. 1994), Salinity measurements from 
the buoys were not available for the location. Thus, an 
isothermal layer depth (ILD) definition is used with a 
AT value of 0.5°C. The ILD can be summarized in its 
simplest form as being the depth at the base of an iso¬ 
thermal layer where the temperature has changed by a 
fixed amount of AT = 0.5°C from the temperature at a 
reference depth of 10 m. A layer depth obtained using 
the ILD with this AT value is approximately equal to 
true mixed MLD in the equatorial ocean (Kara et al. 
2000c). Figure 6 shows the daily buoy MLD as well as 
those from the three NLOM simulations. There are very 
small differences in MLDs from the three experiments. 

Overall the model MLD is shallow. As expected 
based on turbulence considerations (cf. section 3c), the 
MLD in general deepened when subsurface heating was 
permitted within the OGCM. However, the changes at 
this particular location were not significant. For ex¬ 
ample, the ME values between the buoy MLD and those 
obtained from the three experiments are close to each 
other with values of —15.1, -14.8, and -14.8 m, re¬ 
spectively. Large SST differences do not occur during 
winter at (46°N, 131®W) because of the deep MLD (not 
shown). When heat storage is distributed over a thicker 
ocean layer, the amplitude of the SST variations de¬ 
creases (Schneider and Zhu 1998). Since subsurface 
temperatures were unavailable for this buoy, a monthly 
MLD climatology (Kara et al. 2003a) based on tem¬ 
perature and salinity profiles from the World Ocean At¬ 
las 1994 (Levitus et ki. 1994; Levitus and Boyer 1994) 
was used as a proxy for the observed MLD. The MLD 
climatology has been previously shown to be in good 
agreement with hydrographic profiles taken at (49°N, 
13 rW) giving confidence in its use as a proxy for the 
mooring at (46°N, 131°W). 

To investigate NLOM performance in predicting 
MLD we used a normalized rmse defined as 


(nrmse)2 




where v,- (i = 1,2,..., n) is the set of N buoy values 
and Ciit = 1,2 ,..., n) is the set of estimates from the 
NLOM. The nrmse is used because the MLD amplitude 
error tends to be proportional to MLD and use of a 
fractional MLD error prevents the larger errors that oc¬ 
cur when the MLD is deep from dominating the error 
assessment. The nrmse values between the buoy MLD 
and the three experiments are almost identical with val¬ 
ues of 0.33, 0.34, and 0.33. These nrmse values show 
the errors for the NLOM estimates are approximately 
33% of the buoy MLDs. Previously Kantha and Clayson 
(1994) pointed out the influence on the mixed layer of 
upwelling and downwelling processes that are so prev- 
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Fig. 6. Daily averaged MLD time series comparisons between values calculated from buoy 
subsurface temperatures and those from the three NLOM experiments at (0°, 155°W) in 1997. 
See Fig. 5 for monthly values at this location. 


alent at the equator, making MLD simulation from a 
bulk or diffusion model difficult. Taking this into con¬ 
sideration, we can nevertheless state that the NLOM is 
still able to predict MLD well with an annual-mean error 
of only —15 m. 

It is noted that experiments 1 and 2 produce an SST 
that differs only slightly in the annual mean despite the 
strong monthly ^par variation used in experiment 1. This 
very small change in the annual mean SST can be most 
easily explained by examining the compensation depth 


Month 

Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec 




Fig. 7. Compensation depths computed from kp^p for the two moor¬ 
ing locations shown in Fig. 4. Also shown are daily time series of 
model MLD from experiment 1. 


(Dc). The latter is defined as the depth at which the 
PAR decreases to 1%, 


Dc 


In(O.Ol) 

^PAR 


(24) 


and is routinely used in ocean biology (Lalli and Parsons 
1997). In the present context Dc represents the maxi¬ 
mum depth for solar heating of the upper ocean. The 
maximum value of Dc that occurs during the year rep¬ 
resents the full extent of the water column over which 
thermal energy from solar heating is stored. Whenever 
the MLD attains or exceeds the maximum Dc value for 
the seasonal cycle, the thermal energy from solar heating 
that was not absorbed in the mixed layer because of 
penetration below the mixed layer during the remainder 
of the year is entirely entrained into the mixed layer. 
Figure 7 shows the compensation depth computed from 
the monthly values and the MLD from experiment 
1 at the two mooring locations shown in Fig. 4. At both 
locations the model MLD in January deepens to almost, 
or beyond, the maximum value. The same annual in¬ 
tegrated heating is therefore obtained in experiments 1 
and 2 because the mixed layer entrains all of the solar 
heating input into the upper ocean over the seasonal 
cycle. In other words, the mixed layer allows the ex¬ 
periment with space- and time-varying /^p^R to achieve 
nearly the same annual integrated heating as the ex¬ 
periment with constant A:par = 0.06 m"^ Ohlmann et 
al. (1996) previously identified this mechanism when 
they pointed out that penetrating solar energy can be 
trapped within the seasonal pycnocline at mid- and high 
latitudes, with this thermal energy becoming unavailable 
for exchange with the atmosphere until winter mixed 
layer deepening. 
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c. Comparisons with other OGCM studies 

With regard to other studies, we note here that Mur- 
tugudde et al. (2002) report annual-mean SST differ¬ 
ences as large as —0.5°C between their variable and 
constant-attenuation-depth simulations in the Tropics by 
using a primitive equation OGCM that has 19 sigma 
layers beneath a surface mixed layer. They found that 
warming of the water column below the mixed layer 
led to an SST warming in the Tropics despite penetrating 
solar radiation below the mixed layer. In contrast, we 
find the SST decreases when permitting penetration of 
solar radiation. These different outcomes are likely due 
to the different data sources and methods used for cal¬ 
culating the penetration of solar radiation. Murtugudde 
et al. (2002) obtain penetration depths from the CZCS 
data that are from 10--20 m greater than those obtained 
here using the SeaWiFS data. The larger penetration 
depths lead to solar heating being distributed over a 
greater depth of the upper ocean, thereby producing a 
weaker and deeper thermocline. As a consequence, 
wind-deepening events more easily erode the thermo¬ 
cline than when solar penetration is neglected. This 
leads to larger amounts of appreciably warmer water 
being entrained into the mixed layer and, thus, a warmer 
SST and deeper MLD. The lesser penetration depths 
used in this study result in a stronger and shallower 
thermocline. This permits colder water to be entrained 
during wind-deepening events and sustains a shallower 
mixed layer. This combines to produce an SST that is 
colder than when solar penetration is neglected. The 
opposing results from these two studies point to the 
sensitivity of SST in the Tropics to the parameterization 
of the penetration depths. 

Earlier studies consistently find much larger differ¬ 
ences in the mean MLD error between simulations per¬ 
formed when excluding and including sunlight penetra¬ 
tion below the ocean surface than those reported here. 
For example, using a Richardson-number-dependent 
vertical mixing scheme of heat and momentum, Schnei¬ 
der and Zhu (1998) found mean MLD errors increased 
by 30 m off the equator in the Pacific, while the mean 
MLD error increased by 10-30 m at the equator. Mur¬ 
tugudde et al. (2002) obtained annual-mean MLD dif¬ 
ferences of up to 12 m between their variable and con- 
stant-attenuation-depth simulations in the Tropics. This 
is a threefold increase over the annual mean MLD dif¬ 
ference of 2-4 m that we obtained between our corre¬ 
sponding simulations. Nakamoto et al. (2001) report 
mean MLD errors of 20-40 m in the western and central 
equatorial Pacific. These larger differences relative to 
our findings can be related to the greater penetration 
depths obtained using the CZCS versus SeaWiFS data. 
The greater depth penetration obtained from the CZCS 
data permits a larger value for the radiation contribution 
in the ATKE equation [see Eq. (7)] and thus a deeper 
MLD. 

Kantha and Clay son (1994) performed a one-dimen¬ 


sional ocean model study to examine turbidity effects 
on SST. The model was based on second-order closure 
of turbulence. They noted that the SST rise for a very 
shallow MLD is highly sensitive to the parameterization 
of solar extinction. However, their results were limited 
because they used the Secchi depth approach (Gordon 
and Wouters 1978) to determine ocean turbidity in the 
simulations and used it for only one location. They con¬ 
cluded that there could be as large as VC differences 
in SST when using different turbidity in model simu¬ 
lations. This is confirmed by the variety of buoys used 
in this study. Results presented in this paper support the 
conclusion of Kantha and Clay son (1994) that, if avail¬ 
able, a model simulation should include turbidity to ac¬ 
count for variations in solar radiation extinction in the 
mixed layer. However, a constant turbidity (i.e., = 

0.06 m"*) assumption is also adequate to obtain reliable 
SSTs in global ocean simulations. Although there can 
be differences in SST for some months when using a 
variable ocean turbidity rather than using constant ^p^R, 
overall there are no significant differences between SSTs 
obtained from experiment 1 and experiment 2 over an¬ 
nual time scales. This is true for the two buoy locations 
that we examined in detail and statistically for a large 
number of buoys (see Tables 2 and 3). Assuming a 
constant kp^R = 0.06 m“* instead of using a spatially 
and temporally varying /:pAR does not change the SST 
simulation very much. 

4. Summary and conclusions 

In this paper, the impact of solar subsurface heating 
using space- and time-varying turbidity on SST and 
MLD simulations have been examined using the global 
Naval Research Laboratory Layered Ocean Model with 
an embedded bulk-type mixed layer model. The model 
used 6-hourly wind and thermal forcings from the Eu¬ 
ropean Centre for Medium-Range Weather Forecasts 
during 1996 and 1997. However, it develops its own 
sensible and latent heat fluxes from this forcing input 
and model SST. In the model, the attenuation of solar 
radiation with depth was included by using a global 
monthly attenuation of photosynthetically available ra¬ 
diation, denoted as ^par, climatology. These global fields 
of fcpAR are derived from the Sea-Viewing Wide Field- 
of-View Sensor data on the spectral diffuse attenuation 
coefficient at 490 nm This study showed that 

differences in SST and MLD values simulated using a 
spatially and temporally varying attenuation parameter 
versus a constant attenuation value of 0.06 m"^ were 
small when averaged over the course of a year but were 
significant on shorter time scales. 

To determine the impacts of ocean turbidity on SST 
simulations and find an indication of where to use a 
monthly climatology over the global ocean, the 
percentage of PAR absorbed below the mixed layer was 
calculated using global monthly mixed layer depth 
(MLD) and ^par climatologies. In addition, three dif- 
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ferent ocean model experiments were used. The first 
experiment used a spatially and temporally varying 
global monthly itpAR climatology, the second experiment 
used a constant Atrar value of 0.06 m”\ and the last 
experiment assumed all radiation was completely ab¬ 
sorbed within the mixed layer. For the present appli¬ 
cation, the greatest interest is in two situations: 1) where 
subsurface heating occurs below the mixed layer (Fig. 
3), thereby resulting in a cooler SST than obtained by 
assuming complete absorption of solar irradiance within 
the ocean general circulation model (OGCM) mixed lay¬ 
er (expts 1 and 2 vs expt 3) or 2) where sub-mixed 
layer heating is substantially reduced by high turbidity 
(expt 1 vs expt 2). Comparisons of daily SSTs obtained 
from the three NLOM simulations with buoy SSTs 
showed that using monthly varying water type rather 
than using the constant fcp^R approach in all months usu¬ 
ally yields the best SST simulation from the model at 
low latitudes. The SST differences can vary from one 
month to another depending on kp^ and are larger when 
the itpAR value is large. However, considering one year¬ 
long daily SST time series, statistical analysis revealed 
that the differences between experiments 1 and 2 are 
small, at least for the buoy locations examined here. It 
was also found that allowing absorption of all solar 
radiation entering the mixed layer (expt 3) yields de¬ 
graded SST simulations. 

Last, we have also shown the sensitivity of a six-layer 
OGCM to water turbidity and compared it with other 
studies using OGCMs with different vertical resolutions 
and physical parameterizations: a coupled atmosphere- 
ocean general circulation model and a reduced-gravity, 
primitive equation OGCM having 19 sigma layers be¬ 
neath a surface mixed layer. As noted, the effects of 
water turbidity on SST and MLD prediction can vary 
substantially when using different types of OGCMs be¬ 
cause of differences in physical parameterization, forc¬ 
ing, and turbidity product used. 
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APPENDIX 
List of Symbols 


b^w^ 

Buoyancy (m^ s“^) 

^rad 

Contribution due to attenuation of solar ra¬ 
diation (W m"2) 

^pa 

Specific heat of air (1004.5 J kg~^ K"*) 

Gpw 

Specific heat of water (3993 J kg~^ K~^) 

Dc 

Maximum depth of solar heating of the 
upper ocean 

e 

Mean of NLOM SST time series (°C) 

r 

Coriolis parameter at 5° latitude (2.5 X 
10-5 s-i) 

K 

Mixed layer depth (m) 

hi 

Minimum mixed layer depth (10 m) 

hfn 

Mixed layer equilibrium depth (m) 

hp 

Radiation absorption e-folding length scale 
(m) 

^490 

Diffusive attenuation coefficient at 490 nm 

(m-O 

^PAR 

Attenuation coefficient for photosyntheti- 
cally available radiation (m'O 


Coefficient of horizontal temperature dif- 
fusivity (0 m^ 

w,. 

Kraus-Thmer constants (wj = 6.25, = 

7.5, = 6.3, Wg = 0.3) 

ME 

Mean error (°C) 

Wc 

Kraus-Tbmer constant (n^ = 1) 

N 

Number of cases in the time series 

nrmse 

Normalized root-mean-square difference 

P 

Net rate of ATKE generation (m^ s“^) 

PAR 

Photosynthetically available radiation (W 

m-2) 

Qa 

Net heat flux at the ocean surface (W m"^) 

Gir 

Infrared radiation (W m-^) 

Ql 

Latent heat flux (W m-^) 

Glw 

Net longwave radiation at sea surface (W 
m“2) 

Qpiz) 

Penetrating solar radiation at depth z (W 
m-2) 

Qs 

Sensible heat flux (W m"^) 

Gsol 

Solar irradiance at the sea surface (W m"^) 

Qiz) 

Net surface heat flux at depth z (W m~^) 

R 

Linear correlation coefficient 

rmse 

Root-mean-square difference (°C) 

SS 

Skill score 


Time (s) 
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T, 

Temperature just below the mixed layer 
CC) 

T,n 

Sea surface temperature (°C) 


Friction velocity (m 


Layer-1 velocity (m s"*) 

a{TJ 

Coefficient of thermal expansion of sea¬ 
water (°C"0 

An 

Temperature shear (°C) 

Ar? 

Minimum temperature difference (0.2®C) 

An 

Temperature change across the mixed layer 

rc) 


Vertical mixing velocity (m s 0 

a 

Rotation rate of the earth (7.292 X 

s“0 

Latitude (°) 


Po 

Reference density for seawater (1000 kg 

O'e 

Standard deviation of NLOM SST time se¬ 
ries CO 


Standard deviation of buoy SST time series 
CQ 


MLD relaxation ^-folding time (s) 

V 

Mean of buoy SST time series CO 


Wind stress (N m“^) 
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